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ITERATION  METHODS  FOR  NON-LINEAR  PROBLEMS 

ABSTRACT 

The  method  of  successive  displacements  or  relaxation  methods 
are  investigated  for  a  class  of  non-linear  problems.   In  particular 
it  is  shown  that  these  methods  are  applicable  to  a  large  class  of 
non-linear  problems  arising  from  variational  problems  which  yield 
elliptic  equations.   Constructive  existence  and  uniqueness  theorems 
are  presented  for  the  discrete  problem  and  criteria  are  given  for 
a  practical  method  of  obtaining  solutions.   The  example  of  a 
discrete  Plateau  problem  is  used  to  illustrate  the  feasibility 
of  the  results.   The  processes  are  also  shown  to  apply  to  uniformly 
elliptic  problems. 
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ITERATION  METHODS  FOR  NON- LINEAR  PROBLEMS 

1 .  Introduction. 

The  method  of  relaxation  or  the  method  of  successive  dis- 
placements has  been  used  extensively  to  solve  linear  problems  [1]. 
In  particular  this  method  has  been  applied  to  the  solution  of 
linear  elliptic  equations  and  has  been  especially  successful 
with  the  use  of  automatic  digital  computers. 

Analogous  methods  have  been  used  in  practice,  with  apparent 
success,  on  non-linear  problems  as  well.   For  the  most  part, 
these  have  not  been  justified  mathematically  and  this  work  is 
an  attempt  to  fill  this  gap.   In  particular  it  is  shown  that 
the  relaxation  methods  yield  solutions  to  problems  arising  from 
the  minimization  of  certain  convex  functions.   In  practice,  these 
functions  are  obtained  by  approximating  multiple  Integrals  in  a 
calculus  of  variations  problem. 

It  is  shown  that  an  approximate  Plateau  problem  may  be 
solved  by  a  successive  displacements  method,  or  a  method  analogous 
to  Liebmann's  method. 

We  at  the  same  time  obtain  an  extension  of  a  free  steering 
theorem  for  positive  definite  symmetric  matrices  given  as 
theorem  4  of  [3],  and  results  of  Ostrowskl  [2]. 

2.  Definitions. 

Let   E   denote  Euclidean  n-space  and  let   u  be  a  column 
vector  in   E  .   Let   G(u)   be  a  real  valued,  twice  continuously 
dlff erentiable  function  defined  on  a  convex  domain   K  in   E  . 
For  each  u  e  K  we  define  the  matrix   A  =  A(u)  =  (a. .)   by 

u 
O 

a.  .  =  a.  .  (u)  =  ^    G/()u.<)u  .=G        i,j=l,2,...,n. 


-  4  - 


We  assume  further  that  for  all  u  e  K 

a)  A(u)   is  positive  definite 

b)  G(u)  >  m 

m 

for  some  constant  m.   If  u,  v  e  K  then  we  write  w  e  (u,v) 

if  w  lies  on  the  line  segment  joining  u  and  v,  that  is 

if  there  exists  a  number  a,  0  <  a  <  1   such  that  w  =  au  +  (l-a)v. 

Let  r(u)  =  grad  G(u),   then 

(2.1)  G(u)  -  G(v)  =  r(v)'^(u-v)  +  ■|(u-v)^A(w)  (u-v) ,  w  e  (u,v) 

from  the  mean  value  theorem.   (The  superscript   T  denotes 
transpose. ) 

It  then  follows  from  condition  a),  as  is  well  known,  that 
G(u)   is  strictly  convex.   That  is  if  u,  v  e  K  and   0  <  a  <  1 
then 

(2.2)  G(au  +  (l-a)v)  <  aG(u)  +  (l-a)G(v) 

with  equality  if  and  only  if   u  =  v. 
We  seek  solutions  to  the  equation 

(2.3)  r(u)  =  0 

or  stationary  points  for   G(u).   From  (2.1)  and  condition  a) 
it  follows  that  the  solution  is  unique  and  yields  a  global 
minimum  for   G(u).   For  if   u   and   v   are  distinct  solutions  then 

G(u)  -  G(v)  =  ^(u-v)^G(w)(u-v)  >  0 
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and  by  the  same  token   G(v)  >  G(u). 

We  proceed  now  to  give  a  constructive  existence  theorem 
for  (2.3). 

Let   D  be  a  convex  set  such  that   D  C  K,   (the  bar 
denoting  closure),   then  we  call   D  a  solvent  set  if  for 
every  u   e  D  there  exists,  for  each   1,   i  =  l,2,...,n 
a  vector  u'  e  D  such  that 

(2.4)        r.(u°,...,  ul_^,    ul,  u°^^,...,  u°)  =  0 


^k  =  ^k-  ^  ^  ^- 


For  any  given  u   £  D  we  now  define  a  sequence   iu^  i  in 
D  as  follows: 

Let   Z  =  (l,2,...,n)   and  \i    i  be  an  infinite  sequence 

of  integers  chosen  from   Z   such  that  every  integer  in   Z  occurs 

infinitely  often  in  the  sequence.   If  u^  e  D  then  u^  is 
obtained  from  u^   by  changing  one  component  as  follows: 

Choose   Uj\^   '   to  be  the  solution  to   the  equation  in  u. 


(2.5)    r^(u(P),...,  u[P|,  u^,  ^k+1---'  ^n^^)  =  0,  k  =  i 


"P 


while 


u(P+l)  =uiP)  ,   k/  i  , 
k        k    '    '^  p 


The  process  defined  here  may  be  regarded  as  an  extension 
of  the  method  of  successive  displacements  as  used  for  linear 
problems.   We  continue  to  call  the  above  process  by  this  name 
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or  briefly,  a  relaxation  process. 

5.   Convergence  of  the  Sequence. 

We  now  prove  the  following: 

Theorem  3.1.   Let   G(u)   satisfy  the  conditions  stated 
above  and  assume  that  a  solvent  set   D  exists  which  is  bounded. 
Then  (2.3)  has  a  unique  solution  u*   and  the  sequence  \y^  j     of 
the  relaxation  process  converges  to  u*. 

Proof;   The  uniqueness  has  already  been  shown  above. 

We  assume  now  that  x  is  a  limit  point  such  that   r(x)  /  0. 
We  note  from  (2.1)  that 

(3.1)  G(uP)  -  G(uP+l)  =i(uP-uP+l)TA(v)(uP-uP+M,  V  e  (u^^u^^^) 

-  2^^k   ^k     ''  ^kk^'''''       P 
so  that 

(3.2)  G(uP)  >  G(uP+^). 

It  follows  from  condition  b)  that   G(uP)  converges  monotonically 

to  some  number   G_  =  G(x). 

00 

Let   p  =  min.  f|rj^(x)  I  r^{x)   /  0,  i  g  Z^   and  since   D   is 
bounded  there  exist  positive  numbers   A,y\  such  that 


(3.3)  AW  w  <  w  A(u)w  <  Aw  w 


for  all  u  e  D. 

Let   U  be  a  neighborhood  of  x  given  by 
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O-^)  |u-x|  <  6  <  1^ 


where  the  distance  is  Euclidean.   We  choose  an  Integer   N 
sufficiently  large  that  for  all   p  >  N 

(5.5)  ■        G(uP)  -  G(uP+l)  <  ^   . 

Then  for  all  p  >  N  we  get  from  (3-1),  (3-5)  and  (3-5)  that 

G(uP)  -  G(uP+l)  >|  |uP-uP+l|2 

(3.6)  |uP-uP-^l|  <4. 

If  now  for  some  p  >  N,   there  is  a  uP  e  U  such  that 

r.  (x)  /  0   then,  we  get  that,  for  w  e  (uP,x)   and  k  =  i 
P  ^ 

|r^(uP)  -  r^(x)|  =  1^  a^.(w)(Uj(P)-x.)|  <  |a(w)(uP-x)|  <  P/2 

J=l 

and  thus 


k 


(uP)|  <  |r^(x)|  -  Ir,  ruP)|  <  P/2 


or  P/2  <  |r^(uP)|. 


This  implies  that   uP"*"  /  uP   and  that 


P/2  <  |r^(uP)|  =  Ir^(uP)  -  r^(uP+l)|  =  a^^  (v)  |uP-uP+l  1  <^|uP-uP+^| 


which  contradicts  (3.6). 
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We  may  therefore  assume  that  for  all  u^  e  U,   such  that 


>  N,  r^  (x)  =  0  and   (uP'^^-uP)'^r(x)  =  0. 


P 

Let    K  be  the  set  of   u  e  K  such  that   G(u)  <  G(x) 

00  ^  '  _   ^  ' 

and  let   P  be  the  bounding  surface  G(u)  =  G_,  =  G(x)   of 
'  oo  ^  '     oo 

K   .   Let   T  be  the  tangent  hyperplane   (u-x)  r(x)  =  0   to 

P   at  X.  Denote  by  E''"(T)   and  E~(T)   the  open  half  spaces 
defined  by  T  where  E  (T)   Is  the  set   (u-x)  r(x)  <  0,   and 


E  (t)   contains   K 


00 

The  latter  follows  from  the  convexity  of   K    and  (2.1). 
This  also  implies  that   T  has  only  the  point   x  in  common 

with    n  . 

'oo 
Let   T   be  a  hyperplane  parallel  to   T  such  that   T  CI  E~ (t) 

and  is  a  distance   x  >  0  from  T.   Let 


c   =  e'^(t  )  n  K 


00 


then  it  follows  from  the  strict  convexity  of   G(u)   that,  for 

sufficiently  small  i,      C   is  bounded. 

Let   T^  be  so  chosen  and  let   p.   be  a  positive  constant 

such  that   w'^a(u)w  >  |j.w  w  for  all   u  e  c^   and  all   w  £  E^. 

If   0  <  T  <  Tq   then  for  u  e  C   0  >  G(u)  -  G(x)  =  (u-x)  r(x)  + 

+  i(u-x)  A(v)(u-x)   for  some   v  e  (u,x).   Thus 


T  >  -  r(x)'^(u-x)  >  ^   ^x|u-x| 


This  implies  that  for  a  sufficiently  small  neighborhood   U 

of   X,   satisfying  (3-4),   C  CZXi,      for  some   t-,  >  0.   Let 

'^1 
T-,  =  T    and  let   S  =  E  (T-^)  y  U. 
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Since   S  ID  K      ,      we  may  assume  that  all  but  a  finite  number 
of  the  u^   are  in   S.   Let   u^  e  E^(T^)  0  U  then  uP^-""  e   S 
and   (uP^^-uP)'^r(x)  =  0.   Thus   u^^^  e  E'^(T^)   and  therefore 
u^"^   €  u.   Thus  all  but  a  finite  number  of   u^   are  in   U. 
If,  say,   r  (x)  y^  Q.    then   1=1   for  at  most  a  finite  number 
of  values  of  p. 

This  contradicts  the  hypothesis  and  thus   r(x)  =  0.   Since 
the  solution  is  unique  the  sequence  converges,  and  the  proof 
Is  complete. 

4.  Converse. 

A  partial  converse  to  theorem  3.1  can  be  obtained  in  the 
following  sense: 

Theorem  4.1.   Let   G(u)  e  C  ,      assume  that   A(u)   is 
non-singular,  and   aj^^(u)  >  0  for  u  e  Kj  i  e  z.   Let   K  be  an  open 
solvent  set  such  that  for  every  u   £  K   the  resulting  relaxation 
process  converges  to  the  same  solution  u*   of   (2.3)   and  u*  £  K. 

Then   A(u*)   is  positive  definite. 

T 
Proof.   Let   w  A(u*)w  <  0  for  some   w  ^^  0,   then  this  is  still 

true  in  some  neighborhood   N:|u-u*|  <  6.   We  may  assume  without 

loss  that   |w|  <  6   and  let   u°  =  w  +  u*.   Since  u°  e  N,   for 

V  £  (u°,u*),0  >  G(u°)  -  G(u*)  =  |(u°-u*)^A(v)(u°-u*).   But  if 

a  relaxation  process  is  begun  with  u  ,  G(u°)  >  G(u*)   from  (3.I). 

Thus   w'^A(u*)w  >  0  for  all   w  e  e  . 

m 

Since   A(u*)   is  non-singular  w  A(u*)w  =  0  implies  that 
w  =  0. 

Theorems  3.I  and  4.1  may  be  regarded  as  an  extension  of  a 
theorem  of  Reich  (see  [l]  p.  237)  for  the  case  where   A(u)  is 
a  constant  matrix  and  the  ordering   (i  ?  is  cyclic  or  periodic 
in  p. 
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5-   Remarks. 

1.  Let  K  =  E   and  assume  that   K  Is  a  solvent  set.   Assume 
furthermore  that  the  set   Kq  =  (u|g(u)  <  G{v°)j,      for  some   v°  e  K, 
is  bounded.   If  the  existence  of  the  set  D  In  theorem  J>.1   Is 
replaced  by  the  above  hypothesis  this  theorem  remains  valid. 

This  follows  readily  from  (3-1)  which  Implies  that   K^^  may 
be  chosen  for   D. 

2.  Let  K  =  E   and  assume  that  there  exists  a  positive 
constant   A  such  that 


(5.1)  w'^A(u)w  >  Aw'^w 


for  all  u,  w  e  K.   Then  K  Is  a  solvent  set  and  K   Is  bounded 

o 

for  any   v°  e  K.   Thus  by  remark  1.  the  conclusion  of  theorem  J>.1 
Is  valid  with  the  above  assumptions. 

To  show  that   K  Is  a  solvent  set  we  let   u  e  K  and  let, 
for  any   1,   u'  e  K  such  that   u'  =  u°,  k  /  1.   We  seek  to  show 
that  there  exists  a  solution  to  r, (u')  =  0  as  an  equation  In  u. 

Since  a. . (u)  >  A  for  all  u.   It  follows  from 


r^(u')  -  r^(u°)  =  a^^(v)(u^-u°),  v  e  (u°,u' ) 


that  for  u^  >  u"  +  |r^(u")|/A,  r^(u')  >  0  while   r^(u')  <  0 
for  ^i  <  ^1  -  |r-(^°)l/>^-   Thus   r.  (u')   has  a  zero  In   K. 

To  show  that  K^     Is  bounded  we  note  from  (2.1)  that  for 
we  (u,v°),0  >  G(u)  -  G(v°)  =  (u-v°)^r(v°)  +  ^  (u-v°  )^A(w)  (u-v° ) 

or  that   lu-v°|  <  2|r(v°)|/A. 


-  11  - 


6.   Approximate  Relaxation. 

In  the  previously  defined  relaxation  process  It  was  assumed 
only  that  the  equation  (2.4)  had  a  solution.   In  practice  one 
would  like  to  have  a  process  which  Is  more  constructive.   The 
natural  approach  would  be  to  take  several  steps  In  Newton's 
method  to  find  an  approximation  to  the  root.   Since  In  the 
linear  case  only  one  such  step  Is  taken  we  consider  this 
possibility  first. 

Assume  that   K  =  E   and  that,  for   G(u)   satisfying  the  a),  b) 
conditions  of  section  2,    the  set   K  =  (u|g(u)  <  G(u°)j  Is 
bounded,  for  some   u   t  K.   Let 


a.  =   mln   a. . (u) ,  1  £  Z, 

and  assume  that  there  exist  positive  constants   p.   such  that 
a. . (u)  <  p.   for  all   u  e  K. 

Given  a  sequence  o^  Indices    \}-rrS      which  exhaust   Z 
infinitely  often,  and  a  sequence  of  numbers   icu  /,  p  =  0,1,2,..., 
then  for  a  given  u   we  define  an  approximate  relaxation  process  by 


(6.1) 


^(p+1)  ^   (P)  +  ^  h  ,   k  =  1 
1        k      P  P        P 


^(p+1)  =  ^{^)  k  ^  i 

^k       ^k         '   ^  ^  ^p 


where  h   =  -  r,  (uP)/a,  ,  (u^).   For  this  process  we  prove 

p      k  ^   •' '  kk  ^   '  ^ 

Theorem  6.1.   For  a  suitable  choice  of  fco  j   the  approximate 

relaxation  process  (6.1)  converges  to  a  unique  solution  of  (2.5). 


-  12 


Proof.  The  proof  is  similar  to  that  for  theorem  J>.1 
Let   Y„  =  rain (a, /p. )   and  let  y     be  chosen  such  that 
0  <  7  <  Y^  <  1 .   If  the   CO   are  such  that 


0  <  7  <  cJOp  <  2Yq  -  T,   p  -  0,1,2,  .  .  . 
then  from  (2.1),  for  v  e  (uP,uP"^^),  k  =  i  , 


=  -pr  (X)  h    (cu  a,  ,  (v)-2a,  ,  (u^)) 
2  p  p  ^  p  kk  ^  '    kk ''  ^ ' 


<-^Vp^kk(^^)/27o- 


Thus  it  follows  that  the  G(uP)  are  again  monotone  non- 
increasing,  all  the  u^  are  in  K^  and  that  r  (u^)  — >  0. 
The  remainder  of  the  proof  follows  in  the  same  manner  as  for 
theorem  5-1- 

Remark  3.   A  converse  similar  to  theorem  4.1  can  be  formulated 

for  the  approximate  relaxation  process  providing  a.,  p.   exist 

and   0  <  0)  <  27^.   In  fact  theorem  4.1  can  be  formulated  for 
—  p  —   '0 

all  processes   7^  j   which  begin  in  an  open  set  and  yield  non- 
increasing  sequences   iG(u^)j  . 

7.   Application  to  Elliptic  Equations. 

We  now  consider  an  application  of  the  preceding  theorems  to 
the  problem  of  solving  certain  non-linear  difference  equations. 
In  particular  we  consider  those  equations  which  arise  from 
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minimizing  problems  which  are  discretizations  of  variational 
problems , 

Let   f2  be  a  given  bounded  domain  in  the   (x,y)-plane   with 
boundary  Q,     and  closure  "fi.   Let   F(x^y,p,q)   be  a  given  function 
defined  for   (x^y)  e  q     and  all   (p^q)   such  that   F   is  of  class 
C    in  p,q   and 


i)     If     I,    If     I,    If     |<m<co 
PP  pq  qq    — 


F  F 

pp   P^  I 

ll)  J  =  I  I    is  positive  definite 

F  F 

pq   qq 


for   (x,y)  e  q     and  all   p,q.   (Subscripts  in  p,q   indicate 
differentiation  In  those  variables.) 

Let  a  square  lattice  of  mesh  width  h  >  0  be  defined  by 
the  nodes   x  =  |j,h,  y  =  vh  ;  [i,  v  =  o,  +  1 ,  +  2^  .  .  .     .   To  each 
node   P,  =  (x,y)   we  associate  four  neighbors   (x+h,y),  (x,y+h), 
(x-h,y),  (x,y-h)   and  denote  these  by   P^.^,  ?^^,    P^  ,  P^^^ 
respectively.   Denote  by   (1   the  set  of  nodes   P, ,  P^^...,  P 
in  Q,     all  of  whose  neighbors  lie  in  Q.      We  choose   h   suf- 
ficiently small  that   a   is  not  null. 

Denote  by   a   the  set  of  neighbors  of  nodes  of  Q,.       which 

are  themselves  not  in   a  .   We  list  these  as   P   ,,...,  P 

Ti  n+1       n+m 

and  call  them  boundary  nodes.   The  union  of   fi,   and  Q,,       is 
denoted  by  Ti,^^ . 

Let   u(P)   be  a  real  valued  function  on  a   then  we  write 
u.  =u(P.),  i  =1,...,  n+m,   and  for  those  neighbors   P. ^ 
of   P^   which  are  in  T^   let   u^^  =  u(P^^).   Let  P-^,    ?^,  .  .  .  ,    P^ 
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be  the  nodes  of   a   such  that   P.,  and   P.„  are  In  fi  , 
so  that   N  >  n.   For  each  such  node   P,   we  associate  the 
pair  of  numbers  p^  =  (u^^-u^)/h,  q^  =  (u^2"^i)A- 

Let   b(p)   be  a  given  function  defined  on  the  boundary 
nodes  and  let  u(P)  =  b(P)   for  P  e  6  .   Let  u  =  (u-,,Up,, 
and 


.u  ) 


(7.1) 


N 


G(u)  =  h^  ^  ^(Pk'Pk'^k^ 

k=l 


It  follows  readily  from  11)  that  G(u)  Is  strictly  convex. 
We  however  will  prove  the  following  stronger  result  for  the 
matrix  A(u)   as  defined  above  In  Section  2; 

Theorem  7-1.   The  matrix  A(u)   Is  positive  definite  for  all  u. 

Proof.   Let  w  be  any  non-zero  vector  In  E   then  we  show  that 
w^A(u)w  >  0.   Let   F^^^  =  ^^^(^1^^?^^^^)  ^   let 


A 


^^^  =  (a!^h  -  (^ )   and  set   V,  =  w'^A^^^w   then  w'^Aw  =  \~ 

^  1^  J  k=l 


V, 


J 
(k) 


We  note  first  that   a^'^^  =  0  except  for  at  most  a  principal  sub- 
matrix  of  order  three  arising  from  the  three  variables   u,  ,  u,  , , 

^  k   kl' 

u,  p.   For  notatlonal  convenience  we  drop  the  k   and  write  these 
variables  as   u  ,  u  ,  u  .   If   P  ,  P^i  '  ^b-p  ^  %   then  this 
submatrix  has  the  form 


(7.2)   ^ 


(k) 


F   +  2F    +  F 
pp      pq     qq 


-F    -  F 

PP   pq 


■F    -  F 

qq   pq 


-F   -  F 

pp   pq 


pp 


pq 


-F   -  F 

qq   pq 


pq 


qq 
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>2.(k) 


where  0^^^    =    {H-^ \~   ).  ._p,  -,  o   and  the  elements  of  j6^^^ 

are  evaluated  at  i?^,p^,c\.^^)  • 

vCk) 

We  note  first  that  every  principal  subdeterminant  of  p        , 

order  not  greater  than  two,  Is  positive.   That  the  p\-      >  0 
follows  readily  from  11)  and  the  fact  that  the  sum  of  the 
elements  of  a  positive  definite  matrix  is  positive.   Every 
principal  subdeterminant  of  order  two  has  the  value  A  =  det  J 
and  is  therefore  also  positive.   Since  the  row  sums  of  0r    ' 
are  zero,  0r    '      has  rank  two  and  is  therefore  positive  semi- 
definite. 


If  at  least  one,  but  not  all,  of  the  P,  ,  P,  ,  ,  P,  „   are 

k'   kl'   k2 


m 


% 


the 


n     0r    '      is  by  the  above  argi^ment  positive  definite. 


Thus  in  all  cases   V,  >  0,   and  if  for  some   w  /  0..  w  A(u)w  =  0 
then  Vj^  =  0  for  k  =  1,2,...,N. 


Let   P,   be  any  point  in  a   and  suppose  that   P 


kl 


or 


P,  2   is  in  h..      Then,  since  0S    '      is  positive  definite. 


w, 


0.   Now  let   P,  ,  P,  ,  ,  P,^   all  belong  to   a  .   Then  0 


(k) 


is  given  by  (7-2)  and  may  be  factored  uniquely,  by  the  LDU 
theorem,  into  the  form 


0 


(k) 


10    0 
1     0 

1 


10 


c    c 
20   21 


^00 

0 

0 

1 

=  10 

=20 

0 

A 

0 

0 

1 

c 

21 

T 
=   U   DU 

0 

0 

0 

0 

.0 

1 



,  , 

, 

_ 

If  we  use  the  same  local  indices  with  the 


w,   as  with  the 
k 


u,  ,   associate   w.,  w^  ,  w^   with  the  points   P,  ,  P,  -,  ,  P,  ^   and 
k  0   12  ^  k   kl    1 


k2 


;et   w  =  (Wq,  w^,  w^)  ,  Uw  =  V,  V  =  (vq,  v.^  ,  v^ )     then 


— T  — 
v^Dv 


w. 


0.   This  yields   v   =  v.,  =0  and  therefore  w   =  Vp, 


-  C^^  Wp,  Wq 


^10  '^l 


^20  ^2' 


That  is   w   =  0   if   w^ 


0. 
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We  may  then  proceed  to  the  triple  of  points  with  vertex  at 

^2'  ^^   ^^   ^^^      ^1  "  ^k2  ^^^   ^^      ^  ]1   °^   ^12  ^  ^h  ^^^^ 
w.  =  Wp  =  0.   If  not  we  may  repeat  the  above  argument  until 

we  arrive  at  such  a  point  and  the  theorem  Is  proved. 

8.   Example  of  a  Solvent  Set. 

2     2 
Let  F(p,q)  =  f  (cd)   where   cd  =  p   +  q   and  assume  that 

f'(a))  >  0  for  all   cu  _>  0.   If  we  denote  cu,  =  "^(P^)   then 

k=l    ^ 


=  2h^[f'  (uj.4)q^4  +  f '(a^i^)?!^  "  ^  '  (o).  )  (p  . +q^  )  ] 


1  "^ 
=   2hd[u.  -  ^  ^\    6,  u,  ] 


4_ 
a^l 


where   6,   =  f '  (o)  )  ,  a  =  1 ,  2  5  S .   =  f '  (cd.  ) ,  a  =  ^,  ^   and 
la       1  la       la 

4 

/ la 

a=l 

Let   b   =  mln  b(p),  b . .  =  max  b(p)   for   P   taken  over  the 

boundary   fi,  ,   and  let   D  be  the  hypercube   t)^^  <  u^  <  bj^, 

1  =  l,...,n.   Then   D  Is  a  solvent  set.   For  If   u  e  D  and 
for  any   1  £  Z   consider 

r,(u|°),u,)  .Sh5(u,  -i^  «l,u(°)l 


a=  1 

where  we  write  u.  ,   for  the  variables  other  than  uj^   -This 

function  of  the  one  variable  u.   has  the  property 

r^(u_[°V,u^)  >  0  for  u^  =  bj^  while   r^(u^7,u^)  <  0  for  u^  =  b^^. 
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Thus  there  Is  a  zero  on  the  Interval   [b  , bj^.]   and   u'  e  D. 
Thus  we  obtain,  in  particular,  that  the  solution  u*   satisfies 
the  inequalities   b^j^  <  u*  <  b„,  i  e  Z.   If  Q        is  connected 
then  the  equality  will  hold  if  and  only  if   u*(P)   is  constant 
on  H^. 

In  the  case  of  minimal  surfaces   f(cD)  =  \/l  +  cu  and  the 
above  results  will  be  valid.   Since  here 


(8.1)  Giuf   >  h^  y—  03^^  >  h^A^u'^u, 

k-1 

where   A   is  the  smallest  eigenvalue  of  the  corresponding 
discrete  Dirichlet  problem,  we  get  that  for  any  u°  £  E 
K„   is  bounded.   Since,  by  an  argument  similar  to  that  given 
above,   E   is  solvent.  Remark  1  applies,  and  any  initial 
guess   u   may  be  used  to  start  the   iteration. 

Since,  for  any   i  £  Z,  u.   enters  in  at  most  three  terms 
of  the  sum  (7-1)  we  get  that 

(8.2)  r.(u)  ^^^^(F^^)  +  F^^^)  +  F^^'^)) 

h^    p      q      p       q 

(8.5)    a..(u)  =^^    l[F(i)  +2f(^)  +f(^)  +f(^5)  +F(i^)] 
^   ^     11^  ^     du^      ^2^  pp      ^pq     ^qq     ^pp      -^qq   ^ 


For   the   case     f(cjD)    =    \/l   + 


to 


1    +   f  (cD.  )^    -    2p^q.         1   +   Qij        1    +  Pi4 


(8.4)         a..(u)    =  ^ i-i  + 


-?  + 


f(a)i)^  f("3i3)^        f(a3.^)^ 
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we  get  that   aj^j^(u)  <  4   for  all   u  e  E^^.   Thus  all  the  conditions 
of  theorem  6.1  are  satisfied  with  p.  =4. 

It  follows  from  (8.4)  that  If  u°   is  any  guess  then 
a^^(u)  >  4h  /G{u°y   -  a^      and  y°  may  be  chosen  to  be  h  /g(u°)-^. 
Thus  we  find  that  approximate  relaxation  can  be  used  to  solve  an 
approximate  Plateau  problem. 

9.   Uniforirly  Elliptic  Problems. 

Assume  now  that  in  addition  to  the  previous  conditions  on 
F(P>P>q)  we  assume  that  there  exists  a  positive  constant  |j,  such 
that  for  all   (x,y)  e  Ti   ,      and  all   (p,q) 

(9-1)  A(5)  >  ^i  >  0 

where   A  (J)   denotes  the  smallest  eigenvalue  of  J  . 

From  (8.3)  it  follows  immediately  that   a.  .(u)  >  4i_l  and, 
by  the  argument  of   Section  5j  E   is  a  solvent  set.   Condition 
(9.1)  also  implies  that  for  every   v°  ^  E  ,   K^   is  bounded. 
We  will  prove  the  following  stronger  results 

Theorem  9.1.   If   (9-1)   holds  then  there  exists  a  constant 
a  such  that 


(9-2)  A(A(u))  >  a  >  0 


for  all   u  e  E  . 


Proof.   Let   v  e  E^^     be  defined  by   v  =  (p^  jQ-^  Jp2^q2' "  '  " '^N'^N^ 

and  let   H(v)   be  the  right  hand  side  of  (T-l)-   If   B  is  the 

n2    2N 
matrix  of  order  2N  given  by   B  =  (  ^  ? ,  )  ■  .  -,    then  it  is 

easily  seen  that   B  is  the  quasidiagonal  matrix 


T 
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B  =  dlag[J,  ,Jp ,  .  .  .  ,J-j]  ,   where  ^   is  the  second  order  matrix 
corresponding  to   (Pj^^q,^)-   Thus   X(B)  >^  [i  and,  for  any  pali? 
V,  v°  e  E^^,      we  get 

(9.3)     H(v)  =  H(v°)  >  (grad  H(v°) )^(v-v° )  +  ^  |v  -  v° | ^ 

from  the  mean  value  theorem  (2.1). 

Given  a  function  u(P)   on  f2,  ,   which  takes  the  boundary- 
values   b(P)   on  Q,,      we  may  write   v  =  Du  +  c   where   D  is 
a  constant  matrix  and   c   is  a  constant  vector.   Since   |Du| 
is  a  discrete  Dirichlet  norm   it  is  well  known  that  there 
exists  a  positive  constant   p   such  that 


(9.4)  |Du|2  >  p|u|2 


Now  let   u(P),  u  (P)   be  any  two  functions  on  Q,,       and 


'h 


v,  V   the  corresponding  difference  vectors.   Since 


G(u)  =  H(Du  +  c),   the  linear  terms  in  the  respective  expansions 
about  u   and   v   are  equal.   Thus  it  follows  from  (9-3)^ 
(2.1)  and  (9-4)  that  for  some   z  e  (u,  u°) 


(9.5)       (u-u°)Ta(z)(u-u°)  >  ^l|v-v°|2  >  np|u-u°|2 


If  w   is  any  vector  in   E   such  that   |w|  =1   then  choose 
a  neighborhood   U;|u-u°|  <  5   of   u   such  that 


(9.6)  |w^(A(u)  -  A(u°))w|  <  np/2 
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Let   u  =  u°  +  S-  w   then  by  (9-5)   there  Is  a  z  e  U  such  that 
w'^A(z)w  >  M.p.   This  and  (9-6)  Imply  that  w'^A(u°)w  >  ^,p/2 
which  completes  the  proof. 

We  note  from  the  proof  that,  since  the  positive  definiteness 
of   A  was  not  used,  this  yields  an  independent  proof  of  this 
fact  under  condition  (9-1  )• 

That  K   is  bounded  for  every  v   then  follows  from 
Remark  2.   If  we  assume  that  the  largest  eigenvalue   A (5) 
satisfies  the  condition  /\{J)    <  ^   for  all   (x,y)  e  H  and 
all   (p,q)   then  a^(u)  <  4v.   jn  this  case  Theorem  6.1  is 
valid  where  we  may  choose  y     =   m-/^  • 
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